First Principles Calculations of Ionic Vibrational Frequencies in PbMgi/ 3 Nb2/303 
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Lattice dynamics for several ordered supercells with composition PbMg!/3Nb2/303 (PMN) were 
calculated with first-principles frozen phonon methods. Nominal symmetries of the supercells stud- 
ied are reduced by lattice instabilities. Lattice modes corresponding to these instabilities, equilib- 
rium ionic positions, and infrared (IR) reflectivity spectra are reported. 

PACS numbers: 

There is ample evidence of 1:1 (NaCl-type; the "random layer model'—) short-range order (SRO) in PMN2, but 
first principles (FP) calculations with sufficiently large supercells to realistically approximate a SRO-disordered PMN 
crystal are prohibitively time consuming. Relatively small supercells that might reasonably approximation include 
the [001] wee structure which was predicted to be the PMN cation-ordering ground state (GS)£. 

One objective of this study is to compute the relaxed GS for different small supercells consisting of 15 and 30 
ions, to compare their energies, to compute the relaxed coordinates, dynamical charges, vibrational frequencies, and 
infrared (IR) reflection spectra for the GS. These results are important to understand possible contribution of the 
ordered structures to IR reflection and Raman spectra. 

A second objective is understanding the nature of the soft vibrational modes in the ordered structures of PMN. It 
has been reported that relaxors exhibit both ferroelectric (FE) and antiferroelectric (AFE) characteristics, 4 and that 
competition between FE and AFE fluctuations are the cause of glass-type properties in PMN. 

All FP calculations were done with the Vienna ab initio simulation package (VASP)A. Several supercells of PMN 
composition were considered (Fig. flj . Our FP computations show that all these structures are dynamically unstable 
when the ions are placed on ideal perovskite positions, and full relaxation often leads to suprisingly low symmetry. 
For example, the FP GS of the [OOlhviVAf structure (an [001]2:i superlattice) is monoclinic with the Pb close to Mg 
displaced in the (0.18A,0.05A,-0.05A) and (0.18A,0.05A,0.05A) directions. The Pb ions which are between the Nb 
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FIG. 1: PMN structures considered in the present study: [001]jvjva/ (a), [110]jvjvAf (b), [111]jvjvm (c), [IIIJjvt (d), 
[001] NC c>> (e), and [001] (f) 



2 



0.8 
0.6 
0.4 
0.2 
DC 0.0 
0.6 
0.4 
0.2 
0.0 

300 600 900 1200 300 600 900 1200 
v [cm 1 ] v [cm 1 ] 

FIG. 2: Computed reflection spectra for different ordered structures of PMN 
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planes are mostly displaced in the x direction by 0.27 A). 



TABLE I: Ionic coordinates (in A) in the 
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Tetragonal [OOlJjvcc PMN is also dynamically unstable, and has has a wide spectrum of instabilities that are 
associated with FE, AFE, octahedral tilting and other modes. We relaxed the [001] wee structure after a random 
initial perturbation of the ions. Relaxed coordinates arc listed in Table [TJ The basis vectors are: (4.00, -4.00, 0.00), 
(4.01,4.01,0.01), (0.00,-0.01,3-4.07) relative to cubic perovskite. Relaxations of ionic coordinates are a complex mixture 
of octahedral deformations around Mg-ions (with rather large frozen angles) plus AFE, and FE, Pb-displacements in 
opposition to neighboring O-ions. The x-direction displacements differ from y-direction displacements. Pb-ions have 
cither two, or four, Mg nearest neighbors (nn): Pb-ions with two Mg nn are displaced, by 0.25 A, towards the centroid 
of the nearest Mg-Mg second-nn pair. Pb-ions with four Mg-nn are not displaced. Although we can not prove that 
the structure in Table [I] is the PMN [001] ncc" GS (further relaxations may still reduce its energy in the case), it is 
dynamically stable and has lower energy than PMN [001]jvcc that is relaxed with only FE displacements, or only 
octahedral tilting. 

Dynamical charges Z* for the tetragonal [001] jvcc structure are listed in Table ITU Nb ions that have Nb nn 
in the ±a directions have particularly large Z* aa . Dynamical charges for Pb also exhibit significant anisotropy, and 
environment dependence. 
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TABLE II: The dynamical charges of ions in the 30-ion tetragonal supercell of PMN 
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According to Ghosea&, dynamical (Born) charges Z* are associated with the average field E av but not with the 
local field Ei, as the Szigeti charge Z isi However the product of the Born charge with the average field Z*E av must 
equal the product of the Szigeti charge with the local field, ZE\ — . It follows that 

Z* = (1) 
Consider the average field in a stack of xy-layers with different dielectric permittivities £\ and £2: 

E„=DL(± + ±) (2) 

where D is an external field, L is sample length parallel to the field, d\ and d,2 are the lengths of stacked layers. 
The local field in a uniform dielectric is 



Ei = - (3) 

£% 



Substituting these expressions into the initial one yields: 



* Zi 
Z* = j; 



£1 £2 



(4) 



It follows that, with two layers, one with low and the other with high electronic dielectric permittivity, Z* in the 
more polarizable layer decreases, and Z* in the less polarizable layer increases. The FP calculations follow this trend: 
Z* Mg , which is less polarizable, increases while Z*^ b , which is more polarizable, decreases. Within the Nb plane, one 
calculates more normal values for Z^ b . 

For the 30-ion [001]atcc7' structure, computed diagonal frequencies of the dynamical matrix are listed in Table ITTT1 
The frequencies of Nb-, Mg- and bending O-vibrations are all between (200 to 350) cm -1 . The O-vibrations along 
Mg-0 and Nb-0 bonds range from (520 to 735) cm -1 . The Pb diagonal frequencies are in the range (67 to 108) cm -1 . 

As in experiment 8 , computed IR reflection spectra (Fig. |2J consist of three main reststrahlen bands of the vibrational 
modes that are typical of perovskites: the first group is below 120 cm^ 1 ; the second spreads from 150 cm -1 to 400 
cm -1 ; and the third is from 500 cm -1 to 800 cm -1 . The two lower bands split into two subbands each. Assignments 
of these bands can be made on the basis of the diagonal frequencies shown in Table IIIII In experimental dataS, 
these groups of lines are rather compact as in the computed [OOIJatcc and [111] maw structures. However the 
experimentally determined magnitude of the reflectivity in the interval from 500 cm -1 to 800 cm -1 is lower than 
in the computation. This could be connected with an overestimation of the oxygen dynamical charge and/or with 
large damping for some frequencies in this interval (an estimated damping constant of 60 cm^ 1 was used for all 
frequencies), or it could be that the systems studied here are not sufficiently representative of SRO-disordered PMN. 

The lowest calculated optical frequency in the equilibrium 30-ion [001]jvcc structure is 24 cm -1 . It is lower than 
the lowest Pb diagonal frequency (60 cm -1 ) shown in Table ITTT1 because of the interaction among the ionic vibrations. 
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TABLE III: The diagonal frequencies of the dynamical matrix for the 30-ion [001]jvcc supercell of PMN 
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This mode is a mixture of opposite but not fully compensated Pb-displacements. Other Pb-related modes are spread 
over the interval (24 to 129) cm -1 . Ferroelectric Pb displacements are at (40 and 60 to 90) cm -1 although some 
contribution to e exists in the whole interval from 24 cm -1 to 129 cm -1 . Some modes in this interval have large Pb 
components for definite sites that implies that these modes are quasilocal. The same can be said of some modes in 
the frequency interval from 500 cm" 1 to 800 cm -1 : some of these modes have a very large components corresponding 
to oxygen vibrations along Nb-O-Mg bonds. 

Raman spectra show broad lines with gradual temperature dependence in a wide temperature interval2*ifi. The 
presence of these lines would be forbidden if the ions were in symmetric environments. Ionic displacements due to 
disorder and symmetry breaking can explain the existence of these Raman lines. A possible measure of the intensities 
of these lines is the square of the projection of the ionic displacements, from symmetric positions, onto the vibrational 
modes: S% — Kv^Vi)! 2 , where is the vector of the frozen displacements, and Vj is the vector of the i-th vibration 
in the displacements' representation. 

Low-frequency Raman lines (about 50 cm -1 may be related to the quasilocal Pb-0 vibrations 10 . Note that as 
is typical of Pb-based perovskites that also have large B-cations, such as Mg, that Pb-vibrational branches have 
relatively small dispersion and are unstable across most of the Brillouin zone. The particular instabilities that freeze 
in depend sensitively on the local electric fields produced by the Nb-Sc configuration. Freezing of lattice instabilities 
creates low-symmetry Pb-sites, and also allows AFE Pb-vibrations to couple with FE vibrations. 

Calculated structural instabilities in the various supercells imply reductions in their formation energies relative to 
the values reported in Burton and Cockayne^, but the hierarchy of formation energies and predicted GS, [001]atcc"j 
are the same. Column one in Table lTvl corresponds to the structures with symmetry restrictions imposed (for instance, 
tetrgonal symmetry was imposed for [001] m n n and [001] jvcc)- Column two includes relaxation energy, to the GS 
for each structure. 

There are many possible ordered derivatives of the random layer modei that have (111) Nb-layers which alternate 
with (111) (Nbi/3Mg 2 /3)-layers in which Nb and Mg are ordered; e.g. Fig. QJi. Some mixed layers in this structure 
are ordered in stripes and others in an Nb- honeycomb pattern with Mg's at hexagon centers. This structure has a 
large formation energy relative to the others in Figs. Apparently intralayer cation ordering increases the total 
energy relative to a disordered configuration. This is consistent with the conclusion of Hoatson et aliii who obtained 
an improved inverse Monte Carlo fit to NMR data with the assumption of Mg-Mg and Nb-Nb clustering within mixed 
layers. The structure shown in Fig. ^ also has higher energy than [001]jvcc- 

TABLE IV: The supercell's energy (in eV per 30-ion supercell) 
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Finally, our FP computations have shown that the relaxation of the ionic coordinated in the small supercells of 
PMN does not change the hierarchy of the energy of these structures: [001]atcc" remains to be the GS among the 
small supercells considered. The low frequency vibrations in this structure is connected with AFE Pb vibrations and 
with Pb-0 stretching mode. The computed IR reflection spectrum, qualitatively, corresponds to experimental data 8 



although there are some discrepancies in the line magnitudes at high frequencies. 
S.A.P. appreciates discussions with Toulouse, Svitelskiy, Petzelt, Kamba and Yuzuk. 
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